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We consider Higgs boson production through gluon-gluon fusion in hadron collisions. We present 
a numerical program that computes the cross section up to NNLO in QCD perturbation theory. 
The program includes the decay modes H — > yy, H — > WW — > Ivlv, H — > ZZ — > 4 leptons, and 
allows the user to apply arbitrary cuts on the momenta of the partons and of the photons or leptons 
that are produced in the final state. 
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1. INTRODUCTION 



Gluon-gluon fusion is the main production channel of the Standard Model Higgs boson at 
the LHC. At leading order (LO) in QCD perturbation theory, the cross section is proportional 
to of, as being the QCD coupling. The QCD radiative corrections to the total cross section are 
known at the next-to-leading order (NLO) [jl], ||, ||] and at the next-to-next-to-leading order (NNLO) 
[J|, H| ||, 0, ||, ■ The effects of a jet veto on the total cross section has been studied up to NNLO [|l(]] . 
We recall that all the results at NNLO have been obtained by using the large-M ? approximation, M t 
being the mass of the top quark. 

These NNLO calculations can be supplemented with soft-gluon resummed calculations at 
next-to-next-to-leading logarithmic (NNLL) accuracy either to improve the quantitative accuracy 
of the perturbative predictions (as in the case of the total cross section [ )TT| , [12[ ]) or to provide reli- 
able predictions in phase-space regions where fixed-order calculations are known to fail (as in the 
case of the qj distribution of the Higgs boson [13, 14, 15] at small qj). 

A common feature of these NNLO and NNLL calculations is that they are fully inclusive 
over the produced final state (in particular, over final-state QCD radiation). Therefore they refer to 
situations where the experimental cuts are either ignored (as in the case of the total cross section) 
or taken into account only in simplified cases (as in the case of the jet vetoed cross section). The 
impact of higher-order corrections may be strongly dependent on the details of the applied cuts and 
also the shape of various distributions is typically affected by these details. 

The first NNLO calculation that fully takes into account experimental cuts was reported in 
Ref. [16], considering the decay mode H — > yy (see also [17]). In Ref. Jl8| ] the calculation is 
extended to the decay mode H — > WW — > Ivlv (see also [19]). 

In Ref. [20] we have proposed a method to perform NNLO calculations and we have applied 
it to perform an independent computation of the Higgs production cross section. The method is 
completely different from that used in Refs. [16, 18]. Our calculation is implemented in a fully- 
exclusive parton level event generator. This feature makes it particularly suitable for practical 
applications to the computation of distributions in the form of bin histograms. Our numerical 
program can be downloaded from [21]. The decay modes that are currently implemented are H — > 



yy [gCp, H -» WW -» Ivlv wdH^ZZ^A leptons [gg]. 

In the following we present a brief selection of results that can be obtained by our program. 
We consider Higgs boson production at the LHC and use the MRST2004 parton distributions 
with parton densities and as evaluated at each corresponding order (i.e., we use (n+ l)-loop as 
at N"LO, with n = 0,1,2). The renormalization and factorization scales are fixed to the value 
Mr = = M H , where M H is the mass of the Higgs boson. 



2. RESULTS FOR THE DECAY MODE H -> yy 

We consider the production of a Higgs boson of mass M# = 125 GeV in the H — > yy decay 
mode and follow Ref. [^4|] to apply cuts on the photons. For each event, we classify the pho- 
ton transverse momenta according to their minimum and maximum value, pTmm and /?7- max - The 
photons are required to be in the central rapidity region, |tj| < 2.5, with PTmm > 35 GeV and 
Prmax > 40 GeV. We also require the photons to be isolated: the hadronic (partonic) transverse 
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Figure 1: Distributions in prmm and prmax for the diphoton signal at the LHC. The cross section is divided 
by the branching ratio in two photons. 



energy in a cone of radius R = 0.3 along the photon direction has to be smaller than 6 GeV. Using 
these cuts the impact of the NNLO corrections on the NLO total cross section is reduced from 19% 
to 11%. 

In Fig. |] we plot the distributions in prmm and prmax of the signal process gg — > H — > 77. 
We note that the shape of these distributions sizeably differs when going from LO to NLO and 
to NNLO. The origin of these perturbative instabilities is well known [25]. Since the LO spectra 
are kinematically bounded by pj < Mh/2, each higher-order perturbative contribution produces 
(integrable) logarithmic singularities in the vicinity of that boundary. More detailed studies are 
necessary to assess the theoretical uncertainties of these fixed-order results and the relevance of 
all-order resummed calculations. 

In Fig. || we consider the (normalized) distribution in the variable cos d* , where d* is the polar 
angle of one of the photons in the rest frame of the Higgs boson l . At small values of cos0* the 
distribution is quite stable with respect to higher order QCD corrections. We also note that the LO 
distribution vanishes beyond the value cos(9^ ax < 1. The upper bound cos0^ ax is due to the fact 
that the photons are required to have a minimum pj of 35 GeV. As in the case of Fig. [IJ, in the 
vicinity of this LO kinematical boundary there is an instability of the perturbative results beyond 
LO. 



3. RESULTS FOR THE DECAY MODE H -> Ivlv 

We now consider the production of a Higgs boson with mass Mh = 165 GeV in the decay 



mode H — > Ivlv. We apply a set of preselection cuts taken from the study of Ref. [g6j. The charged 
leptons have pj larger than 20 GeV, and |t]| < 2. The missing pj is larger than 20 GeV and the 
invariant mass of the charged leptons is smaller than 80 GeV. Finally, the azimuthal separation of 
the charged leptons in the transverse plane (A0) is smaller than 135°. By applying these cuts the 
impact of the NNLO corrections on the NLO result does not change and is of about 20%. 

In Fig.|] we plot the A0 distribution at LO, NLO and NNLO. As is well known [^J, the 
charged leptons from the Higgs boson signal tend to be close in angle, and thus the distribution is 



We thank Suzanne Gascon and Markus Schumacher for suggesting the use of this variable. 
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Figure 2: Normalized distribution in the variable cos 0*. 
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Figure 3: Normalized A0 distribution at LO, NLO, NNLO. 



peaked at small Atp . We note that the effect of the QCD corrections is to increase the steepness of 
the distribution, from LO to NLO and from NLO to NNLO. 

We finally apply the set of selection cuts designed to isolate the Higgs boson signal [p6|]. 
The transverse momenta of the two charged leptons are classified according to their minimum and 
maximum value, prmin and prmax- They should fulfil prmin > 25 GeV and 35 GeV < pjmax < 50 
GeV. The missing pj of the event should be larger than 20 GeV and the invariant mass of the 
charged leptons should be smaller than 35 GeV. To exploit the steepness of the Atp distribution 
shown in Fig. [| should be smaller than 45°. Finally, there should be no jets with rfj larger 
than a given value pf l °. 

The corresponding cross sections, for different values of p^ to , are reported in Table [T[ 
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p v T eto (GeV) 


Onlo (ft>) 


<?nnlo (fb) 


no veto 


21. 26 ±0.05 


22.21 ±0.32 


40 


18.62 ±0.05 


17.38 ±0.34 


30 


17.18±0.05 


15.74±0.35 


20 


14.42 ±0.05 


11.31 ±0.38 



Table 1: Cross sections for pp — > H +X — > W+X — > ZvZ V +X of f/ie L//C when selection cuts are applied 
for different pj et °. 



We see that, even without the jet veto, the impact of radiative corrections is strongly reduced 
with this choice of cuts. The jet veto further reduces the effect of the NNLO corrections, which 
become negative for p™ to <40 GeV. 

4. SUMMARY 

We have illustrated the results of a calculation of the Higgs boson production cross section 
at the LHC up to NNLO in QCD perturbation theory. The calculation is performed by using 
the numerical program HNNLO, which includes all the relevant decay modes of the Higgs boson, 
namely, H — > 77, H — > WW — > Ivlv and H — > ZZ — > 4 leptons. The program allows the user to 
apply arbitrary cuts on the momenta of the partons and of the photons/leptons produced in the final 
state, and to obtain the required distributions in the form of bin histograms. We have presented 
a brief selection of numerical results for the decay modes H — > 77 and H —* WW — > Ivlv. More 



detailed results for the decay modes H — > WW and H — > ZZ can be found in Ref. [22]. The fortran 



code HNNLO can be downloaded from [21] 
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